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ABSTRACT 

Cross-correlation of the DIRBE 100 fim survey with previously published Ha maps 
tests correlations between far-infrared dust and the warm ionized interstellar medium 
in different regions of the sky. A 10° x 12° patch at Galactic latitude b = —21° shows 
a correlation slope Oq = 0.85 ±0.44 Rayleighs MJy“^ sr signihcant at 97% conhdence. 
A set of Ha images over the north celestial polar cap yields a weaker correlation slope 
Oq = 0.34 ± 0.33 Rayleighs MJy“^ sr. Combined with observations from microwave 
anisotropy experiments, the data show roughly similar correlations on angular scales 
0?7 to 90°. Microwave experiments may observe more emission per unit dust emission 
than are traced by the same structures observed in Ha. 
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1 Introduction 


Observations of anisotropy in the cosmic microwave background are complicated by 
the presence of foreground Galactic emission along all lines of sight. At high latitudes 
(l&l > 10°), diffuse Galactic emission is dominated by optically thin synchrotron, 
dust, and free-free emission. In principle, these components may be distinguished by 
their different spatial morphology and frequency dependence. In practice, there is no 
emission component for which both the frequency dependence and spatial distribution 
are well determined. Synchrotron radiation dominates radio-frequency surveys, 
but the spectral index steepens with frequency and has poorly-determined spatial 
variation (Bennett et al. 1992, Banday & Wolfendale 1991). Dust emission dominates 
far-infrared surveys, but its spectral behavior at longer wavelengths depends on the 
shape, composition, and size distribution of the dust grains, which are poorly known 
(Desert, Boulanger, & Puget 1990). Free-free emission from electron-ion interactions 
has well-determined spectral behavior but lacks an obvious template map: free-free 
emission never dominates the high-latitude radio sky, while other tracers of the warm 
ionized interstellar medium (WIM) such as Ho; emission, pulsar dispersion measure, 
or N II emission are either incomplete, undersampled, or noise-dominated (Bennett 
et al. 1992, Reynolds 1992, Bennett et al. 1994). 

The problem is particularly acute for free-free emission: lacking an accurate 
template for the spatial distribution, estimation of free-free emission requires pixel-by- 
pixel frequency decomposition which significantly reduces the signal to noise ratio of 
the desired cosmological signal. Recently, Kogut et al. 1996a proposed using infrared 
emission from diffuse cirrus as a tracer of the diffuse ionized gas responsible for free- 
free emission. Gross-correlation of the GOBE Differential Microwave Radiometer 
(DMR) maps at 31.5, 53, and 90 GHz with the Diffuse Infrared Background 
Experiment (DIRBE) far-infrared maps at 100, 140, and 240 /xm shows statistically 
signihcant emission in each microwave map whose spatial distribution on angular 
scales above 7° is traced by the far-infrared dust emission (Kogut et al. 1996b). 
The frequency dependence of this emission, rising sharply at long wavelengths, is 
inconsistent with the expected microwave dust emission. A 2-component £t to dust 
plus radio emission, 

( U \ /^dust / ly \ ^radio 

- ) -Bj/(Tdust) + Aradio ( - ) 

with 1.5 < /3dust < 2 yields spectral index /dradio = —O.Ot^J for the unknown 
component, strongly suggestive of free-free emission (/3fj = —0.15). At high latitudes 
(|6| > 20°), spatial fluctuations in the diffuse synchrotron emission are uncorrelated 
with dust emission, leaving free-free emission (thermal bremsstrahlung) from the WIM 
as the only plausible alternative: on large angular scales, at least 30% of microwave 
emission from the WIM is traced by the diffuse infrared cirrus (Kogut et al. 1996b). 
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The detection of correlation between dust and ionized gas at high latitudes is 
consistent with the correlation observed between dust and free-free emission along the 
Galactic plane (Broadbent, Haslam, & Osborne 1989). Two questions of interest for 
cosmological observations are the extent to which infrared dust reliably traces high- 
latitude radio emission, and whether the correlation depends signihcantly on angular 
scale. Ha emission is a widely used tracer of the WIM. Recently, McCullough 1997 
demonstrated a correlation between Ha emission and the IRAS infrared cirrus at 
1° angular scale, with an amplitude comparable to the value derived from DMR at 
much larger angular scales. In this paper, I present results from cross-correlations of 
high-latitude Ha maps with the DIRBE 100 pm intensity on angular scales 0?7 to 
10 °. 

2 Analysis 

If Ha emission is correlated with far-infrared dust, the Ha intensity in each pixel may 
be written as 

Hi = aDi + Ui, ( 1 ) 

where H is the intensity in the Ha map, D is the far-infrared intensity in the dust map, 
n represents an uncorrelated component (noise, systematic artifacts, or uncorrelated 
emission features), and i is a pixel index. The presence of correlation between two 
data sets is often quantihed using the Pearson correlation coefficient. 

While this has the advantage of being a simple, dimensionless quantity that does not 
require detailed knowledge of noise or artifacts in the two data sets, it can not show 
whether the physically more interesting correlation slope a varies between different 
data sets at different locations or angular scales. An alternative approach uses the 
correlation slope a, dehned as 




with statistical uncertainty 


6a 


ij 


- 1/2 




where M is the covariance matrix for emission between pixels i and j in the Ha 
map. In the limit Mjj = where a is the instrumental noise, equation ^ reduces 
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Figure 1: Reynolds 1980 Ha map and DIRBE 100 /im map in same sparse sampling. 
Full scale is 1-12 Rayleighs in Ha and 3.7-10.3 MJy sr“^ for DIRBE, with white 
regions brighter for both maps. The arrow shows the direction of lower Galactic 
latitude. 


to a least-squares £t. Note that for totally uncorrelated maps {HD) = 0 averaged 
over many pixels (or many maps), so that (a) = 0 as expected when no correlation 
is present. Equation has the additional desirable property of a well-dehned 
distribution: given two uncorrelated data sets x and y with unit normal distribution, 
the distribution a/5a will also have zero mean and unit standard deviation. The 
presence of a true correlation slope Oq in the parent population shifts the mean htted 
value (a) = Oq but does not alter the width of the distribution. 

Few measurements of diffuse Ha emission exist at high Galactic latitude. 
Figure shows the Ha intensity within a 10° x 12° patch at (/, h) = (144, —21) 
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Fluctuations in 100 fiin Intensity (MJy/sr) 


Figure 2: Reynolds 1980 Ha map and DIRBE 100 pm map plotted point by point. 
The fitted slope oq = 0.85 Rayleighs MJy“^ sr. 

measured by a Fabry-Perot spectrometer with a 0?8 beam sampled on a 1° grid 
(Reynolds 1980). The 100 pm emission measured by DIRBE at 0?7 angular resolution 
is shown in the same sparse sampling. Figure ^ plots the two maps point by point 
and shows the fitted correlation slope a for uniform pixel weight. Quantihcation 
of the correlation slope between the two data sets suffers from two problems. 
Uncertainties in the subtraction of geocoronal emission can contribute systematic 
artifacts comparable to the instrumental noise in each line of sight (Reynolds 1980). 
The small number of independent pixels (121) creates appreciable sample variance. 
The far-IR dust emission has a steep power spectrum, P oc where i oc 6~^ is 
the spherical harmonic order (Gautier et al. 1992, Kogut et ah 1996a). Gradients in 
the Ha map may thus correlate with large-scale features in the far-IR dust emission 
simply by chance alignment. 

These problems may be addressed using Monte Garlo techniques, repeating 
the £t in equation ^ after replacing the DIRBE patch Jdirbe with a set of “control” 
patches selected from distant portions of the sky. Since the same Ha map is used for all 
trials, the distribution of fitted slopes automatically accounts for instrumental noise, 
systematic artifacts, and chance alignment, and provides a simple, robust description 
of the probability for an uncorrelated patch of cirrus to produce a slope a as large 

^Throughout this paper I use the 100 map for its high signal to noise ratio and to facilitate 
comparison with the IRAS 100 /rm survey at hner angular scales. I obtain similar results using the 
DIRBE 140 /im and 240 /rm surveys. 
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(or larger) than the value from the correct part of the sky. Note that Eq. ^ is only 
asymptotically unbiased: although the Monte Carlo technique correctly accounts for 
the width of the distribution in a, instrumental noise in the “template” DIRBE map 
reduces the mean value by a factor + n^), where S and n are the signal and 

noise in the DIRBE map. For the data in this paper, S'/n S> 1 and the correction is 
negligible (less than 1%). 

Far-IR emission is markedly anisotropic, so I restrict the control patches to 
those lying within the latitude range of the Ha map. For computational ease, I center 
each control patch on the pixel centers of the COBE quadrilateralized sky cube at 
resolution index 4 within the range 15° < |6| < 30°, resulting in 64 independent cross¬ 
correlations. To allow for possible correlated esc(5) emission (with no smaller-scale 
correlations), I determine the distribution of htted slope a with the control patches in 
random orientation with respect to the Ha map, and again with the control patches 
hxed with respect to the Galactic plane. 

Figure shows the distribution of the htted correlation slopes a between 
the Reynolds Ha map and the set of randomly oriented control patches. The 
correlation slope oq between the Ha map and the DIRBE 100 /xm map in the 
correct position is also shown. Only 3% of the uncorrelated control helds had 
a > qq, suggesting a marginal detection of correlated emission between Ha and 100 
pm dust, Oq = 0.85 ± 0.44 Rayleighs MJy“^ sr within the Reynolds held. Fixing 



Correlation Coefficient (R/[MJy/sr]) 


Figure 3: Histogram of htted correlation slope between the Reynolds 1980 Ha map 
and a set of 64 DIRBE “control” helds. The correlation oq between the Ha map and 
the DIRBE intensity along the same lines of sight is also shown. 


6 







the orientation of the control patches with respect to the Galactic plane does not 
appreciably alter the results: the distribution of control patches still has zero mean 
and similar width, demonstrating that the observed correlation is not purely an 
artifact of general esc(6) trends in both data sets. Repeating the analysis with the 
Pearson correlation coefficient p instead of the correlation slope a provides similar, 
slightly more significant, results: p = 0.48 ± 0.22, significant at 98% confidence. 

Gaustad, McGullough, & Van Buren 1996 measured Ha emission at 0?1 
angular resolution using narrow-band GGD images in a series of 7° x 7° fields covering 
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Figure 4: Ha mosaic over north celestial pole (Gaustad et ah 1996) and fluctuations 
in DIRBE 100 pm. intensity. A mean has been removed from each data set. White is 
brighter for both maps. 
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the north celestial polar cap (declination S > 81°). Visible structure within the Ha 
helds is dominated by systematic artifacts (ghosts and stellar residuals), with intrinsic 
anisotropy less than 1.3 Rayleighs (Gaustad, McCullough, & Van Buren 1996). To 
reduce systematic artifacts, I reject all pixels in each held lying more than 3 standard 
deviations from the mean in that held, and then re-bin the remaining data into DIRBE 
pixels. Figure ^ shows a mosaic of the re-binned Ha map, along with the DIRBE 100 
fim data. 

Artifacts associated with held edges are evident in the Ha mosaic, and direct 
held to held comparison in regions where 2 or more helds overlap conhrms that much 
of the large-scale structure is systematic in origin. Since power in the far-IR emission 
also increases on larger angular scales, the presence of large-scale systematics may 
be partially ohset by increased sensitivity to dihuse structure in the DIRBE maps. 
To allow for this possibility, I cross-correlate the Gaustad et al. Ha data against the 
DIRBE data using two techniques. The hrst method simply correlates the full Ha 
mosaic with the DIRBE 100 pm map, using a set of control helds at 20° < |6| < 33° to 
estimate the uncertainty. The full mosaic shows no statistically signihcant correlation: 
28% of the control helds had a correlation slope larger than the actual helds in the 
correct orientation, oq = 0.21 ± 0.43 Rayleighs MJy“^ sr. 

The correlation from the full mosaic is dominated by systematic artifacts at 
the held edges. A second method avoids these artifacts by correlating each of the 
sixteen 7° x 7° helds independently, then using the arithmetic mean of the 16 htted 
correlation slopes a* to estimate the mean correlation over the full polar cap. As 
before, a set of control helds provides an estimate for the uncertainty in the mean 
correlation. 17% of the control helds had mean correlation larger than the actual 
helds in the correct orientation, yielding a correlation oq = 0.34 ± 0.33 Rayleighs 
MJy“^ sr on angular scales 0?7-7° within the north celestial polar cap. The results 
do not change appreciably if the control patches are held in hxed orientation with 
respect to the Galactic plane. 

3 Discussion 

Figure plots the correlation slope between the WIM and 100 pm dust emission 
as a function of angular scale. Since the dust power spectrum varies as 6^, the 
data are shown at the largest angular scale probed by each experiment. Data from 
microwave experiments (DMR:Kogut et al. 1996b; Saskatoomde Oliveira-Gosta et 
al. 1997) are expressed in Rayleighs assuming a spectral index j3s = —0.15 and a 
conversion 2 pK/Rayleigh at 53 GHz appropriate for gas at 8000 K. There is no 
compelling evidence for variation in the correlation slope as a function of angular 
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Figure 5: Correlation slopes a between the warm ionized interstellar medium and 
far-infrared dust emission as a function of largest angle probed by each experiment. 
Open symbols refer to Ha data and hlled symbols refer to microwave measurements. 
DMR: Kogut et ah 1996b; Saskatoon: de Oliveira-Costa et al. 1997; McCullough: 
McCullough 1997; all others:this work. 


scale. 

The limited data in Figure ^ appear to fall in two groups, with correlations 
derived from microwave observations lying higher than direct Ha correlations, even 
in the NCP held where both microwave and Ha data are available. Although the 
uncertainties on any individual measurement are appreciable, the mean slope derived 
from microwave observations lies 2.7 standard deviations above the mean of the three 
direct Ha observations. Both observations detect structure correlated with the spatial 
distribution of the infrared cirrus, but the microwave experiments appear to observe 
more emission per unit dust than are traced in Ha. 

Further evidence for “extra” microwave emission traced by the infrared cirrus 
comes from the MAX5 microwave anisotropy experiment (him et al. 1996). Data from 
the MAX5 experiment show statistically signihcant correlations between the IRAS 100 
/im map and millimeter-wavelength anisotropy measurements centered at at 105, 180, 
270, and 420 GHz. The frequency dependence is consistent with a superposition of 
dust plus a second component identihed as either free-free emission or anisotropy in 
the cosmic microwave background; although free-free emission is preferred, the data 
do not distinguish decisively between these alternatives. The amplitude of the second 
component is large: if it results from free-free emission, the implied correlation slope 
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between the WIM and dust is a = 37±19 Rayleighs MJy“^ sr on angular scales 6* ~ 2° 
within the MAX5 /i Pegasi held (him et ah 1996). Such a result would imply even 
greater variation in WIM/dust correlation than already evident in Figure ^ (since 
the MAX5 result can not be attributed dehnitively to Galactic emission, it is not 
plotted). Variation in the electron temperature within the WIM would change the 
ratio of microwave free-free emission to Ha emission. However, the temperature of 
the high-latitude cirrus (presumably heated by the same interstellar radiation held 
as the WIM) is constant within 10% when averaged over broad regions (Reach et 
ah 1995), making problematic any preferential heating of the WIM. Other possible 
emission mechanisms that would produce long-wavelength emission correlated with 
emission from the IR cirrus include hat-spectrum synchrotron from supernovae shock 
fronts and radio emission from rotating charged dust in the WIM (Ferrara & Dettmar 
1994). 

4 Summary 

Cross-correlation of the DIRBE 100 /xm survey with previously published Ha maps 
yields positive correlation with amplitude comparable to previously published Ha 
and microwave results. The combined data show no dependence on angular scale, 
indicating that existing infrared maps can be used to trace microwave emission over 
scales 0?7 to 90°. Microwave experiments show a marginally signihcant trend toward 
observing more microwave emission per unit dust emission than are traced by the 
same structures observed in Ha. Additional observations, both Ha and microwave, 
are required to determine if this results from a statistical huke, variations within 
the WIM, or the existence of competing emission mechanisms traced by the spatial 
distribution of high-latitude dust. 


Gary Hinshaw provided helpful conversations. This work was supported in 
part by NASA RTOP 399-20-61-01. 
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